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Figure 1. Closed sequence of chemical equilibria between X, Y, Z and V. 


1. Introduction 


Consider a closed chemical system composed of four coexisting chemical species denoted 
by X, Y, Z and V, which represent four possible states of a macromolecule operating in 
a reaction network far from equilibrium. As discussed by Wyman 1 ^, such a reaction 
can be modeled as a “turning wheel” of one-step transitions of the macromolecule, 
which circulate in a closed reaction path involving the four possibles states. The turning 
wheels have been proposed by Di Cera et al. as a generic model for macromolecular 
autocatalytic interactions. 

While Di Cera’s model considers unidirectional first order interactions, Murza et al. 
in [3 consider a closed sequence of chemical equilibria. In their approach the reaction 
rates are defined as functions of the time dependent product concentrations, multiplied 
by their reaction rate constants. This type of reaction rates has been introduced in 
Wyman’s original paper 181. 


Following the closed sequence of chemical equilibria in 1^ , the autocatalytic chemical 
reactions between X, V, Z and V (see Figure[T]) are governed by the following 4-parameter 
family of nonlinear differential equations 


X = x{kiy — k^v), 
y = y{k 2 Z - kix), 
z = z{kzv - k 2 y), 

V = viki^x — k^z). 

Functions x{t),y{t), z{t) and v{t) are concentrations at time t of the chemical species 
X,Y,Z,V respectively. The parameters ki = ki 2 — kn for i = 1,2,3,4 are differences 
of pairs of reaction rate constants corresponding to each chemical equilibrium. It can 
be easily seen that system m is identical to Di Cera’s model restricted to n = 4, see 
equation (7) in [^. In that work, Di Cera claims that this family exhibits self sustained 
and conservative oscillations only when the parameter k = (fci, ^ 2 , k^, k/C) is in the three 
dimensional manifold 5 = {k € \ {0} : fci/cs — ^ 2^4 = O}. 

Assuming that the conservation of mass x + y + z + v = 1 applies to the macromolec- 
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ular system its kinetic behaviour is described by the three-dimensional system of 
polynomial differential equations 

X = x{kiy — k4{l — X — y — z)), 

< y = y{k2Z-kix), ( 2 ) 

_ i = z{-k2y + ksil -x-y- z)), 

restricted to the flow-invariant bounded region T = {x > 0, y > 0, z > 0, x + y + z < 1}. 

The system ([2|) is a particular case of the class of three-dimensional Lotka-Volterra 
systems (TVS) 


Xi = Xi 


CLi + bi 


i = 1,2,3, 


1=1 


which has been extensively studied starting with the pioneer works of Lotka 1^ and 


Volterra IT]]. These systems have multiple applications in biochemistry. For instance, 
enzyme kinetics 18] , circadian clocks 14] and genetic networks [l|, Q often produce 
sustained oscillations modeled with TVS. 

Solutions of TVS cannot, in general, be written in terms of elementary functions. So 
that the search for invariant manifolds, first integrals or/and integrability conditions can 
be useful to the analysis of the flow. This approach has found increasing popularity over 
the last few years after the works of Christopher and Llibre [^-0] j which are based on the 
Darboux’s theory of integrability, see for instance 0, li, 12] and references therein. Unfor¬ 
tunately for systems of dimension greater than 2 the behaviour of the flow is not entirely 
known, even when the system is integrable. Of course in the case of non-integrable TVS 
the lack of knowledge is higher and other tools are required. Some results concerning 
the existence of limit cycles for special parameter sets can be found in 10, lla . Il9l . |2C|. 


Additional results about the number of limit cycles which can appear after perturbation 
are presented in Qj. 

In this paper we deal with the global analysis of the flow of system ([2]) restricted to the 
region T. Note that the boundary dD of the region is a three dimensional simplex which 
is invariant by the flow. This boundary is formed by the union of the following invariant 
subsets: the invariant faces X = {(0, y, z) : y > 0, z > 0, y -|- z < 1}, T = {(x, 0, z) : x > 
0, z > 0, X -|- z < 1}, Z = {(x, y, 0) : X > 0, y > 0, X -|- y < 1} and S = {(x, y, z) : x > 
0, y > 0, z > 0, X -|- y -|- z = 1}; and the invariant edges TZy^ = {(x,0,0) : 0 < x < 1}, 
T^xz = {(0, y, 0) : 0 < y < 1}, TZ^y = {(0,0, z) : 0 < z < 1}, TZp^ = {(0, y, 1 - y) : 0 < y < 
1}, TZpy = {(x, 0,1 — x) : 0 < X < 1}, and TZpz = {(x, 1 — x, 0) : 0 < x < 1}. We remark 
that the edges TZxz and TZpy are closed segments formed by singular points. 

In order to make easier the analysis we consider the following subsets in the parameter 
space: S~ = {k G : kik^ — ^ 2^4 < 0}, S~^ = {k G : kik^ — ^ 2/04 > 0}, AfZ = {k G 
: kik 2 k^k 4 / 0} and T’S = {k G : kik 2 > 0, kik^, > 0, kik^ > 0}. We note that S~ 
and 5"*“ together with S (defined above) form a partition of the parameter space R^. We 
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Figure 2. Behaviour of the flow of system ([2]) for k G VS ("15: (a) in the interior of T, which 
contains periodic orbits and the segment TZ formed by singular points; (b) at the faces which 
exhibit an edge TZpy formed by singular points, and heteroclinic orbits giving rise to heteroclinic 
loops; (c) at the faces X\JZ which exhibit an edge TZxz formed by singular points, and heteroclinic 
orbits giving rise to heteroclinic loops. 


also note that the parameter set VS is a subset of J\fZ. 

The main result of the paper is summarised in the following theorem. 

Theorem 1.1 


(a) Suppose that k G VS IH S. 
(a-1) The open segment 


R = 


f ks k4- {ki + k3)z \ r, ^ ^ ^4 1 

V fc4 k4: + kl J fcs + /C4 i 


is contained in the interior of T and every point in R is a singular point. 

(a-2) Let p be a point contained in the interior ofT but not in R. Then the orbit 7 p 
through the point p is a periodic orbit. 

(a-3) Each of the two limit sets of every orbit in y WE is a singular point contained 

in the edge TZpy. Moreover, given two orbits 71 C T and 72 C S such that 

a;( 7 i) = q;( 72 ), then ^( 72 ) = a(7i). 

(a-4) Each of the two limit sets of every orbit in XU Z is a singular point contained 

in the edge TZxz- Moreover, given two orbits 71 C <T and 72 C .H such that 

a;( 7 i) = 0 ( 72 ), then ^( 72 ) = a(7i)- 

(b) Suppose that k 0 VSnS and k 7 ^ 0. The limit sets of every orbit in T are contained 
in the boundary of T and these limit sets are non-periodic orbits. 

From Theorem ll.lf al it follows that when k G 7^51^5 the behaviour of the trajectories 
in the whole region T is illustrated in Figure [ 2 l 
Theorem ll.ll bi is partially a corollary of a general theorem for Lotka-Volterra systems 
stating that no limit points are in the interior of the non-negative orthant if there is no 
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singular points in the interior. See Theorem 5.2.1 (p. 43) in 11[. 

On the other hand Theorem 11.11 completely characterizes the region in the parameter 
space where the corresponding system ([2|) exhibits self sustained oscillations. Thus the 
necessary conditions k € <S for the existence of such behaviour, given by Di Cera et al. [^, 
are here completed with the necessary and sufficient condition k G VS n5. Furthermore 
this oscillating behaviour in the interior of T extends to a heteroclinic behaviour at the 
boundary. Therefore the period function defined in the interior of T is a non-constant 
function; it grows when approaching the boundary. As a final remark, we note that 
oscillations in system ([2|) take place in a parameter set of measure zero, and they are 
only “conservative” ones; i.e. are not isolated in the set of all periodic orbits. 

In dimension greater than two continuous dynamical systems may present chaotic 
motion, in the sense that the distance of points on trajectories starting close together 
increases on an exponential rate. This is not our case. The dynamic behaviour of family 
m is very simple and non-strange attractors appear. In fact, as shown in Theorem 
ll.If bl in absence of periodic orbits every orbit goes from one side of the boundary of T 
to another. Nevertheless, we can remark certain singular situations related to the form 
and location of the limit sets. One of these limit set configurations is described in the 
next result. Before stating it we consider the following singular points in the edges TZpy 
and TZxz, respectively 


Ppy ~ ( k!+k2 ’ fci%2) ’ ^py ~ ( 


ks n k. 


k3-\-k4 


^3 + ^4 I ’ 


( 3 ) 




0 

fcl+fc4’^^ ’ 


Qx2 - ( 0 


’ /C3H-/C2 


,0 


When k is in the manifold VS fl S, the points Ppy and qpy are equal and they coincide 
with one of the endpoints of the segment R defined in Theorem 1 1.1 f a- lb Similarly, the 
points Pxz and are also equal and they coincide with the other endpoint of R. On the 
other hand, when k G VS \ S, we define the following segments contained in the edges 
TZpy and TZxz, respectively 


^py ~ {Ppy Ppy) • £ [0,1]} , 

Sxz - {PxZ 4“ p{qxz Pxz) ■ P £ [0, 1]} . 


To clarify the exposition of the next result we introduce the subsets VS^ = {k G "P : 
ki > 0} and VS- = {k G P : /cj < 0} which form a partition of VS. 

Theorem 1.2 Suppose that k G VS \ S. 

(a) Each of the two limit sets of every orbit in the interior ofT is formed by a singular 
point contained in the segments Spy and Sxz- In particular, given a point p in the 
interior ofT, i/k G P5+ n<S+ or k G VS- n<S“, then a( 7 p) G Sxz and w( 7 p) G Spy] 
and i/ k G P<S+ O S~ or k G VS- D <S+ then a(7p) G Spy and w( 7 p) G s^z- 

(b) Each of the two limit sets of every orbit in yuT, is a singular point contained in the 
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Figure 3. Behaviour of the flow of system ([2]) for k G VS^ fl 5+ or for k G VS^ fl S~ : (a) 
the limit set of every orbit in the interior of T is a point contained in the segments s^z and Spy, 
respectively; (b) at the faces 3^ U S where the orbits are heteroclinic orbits connecting singular 
points located at the edge TZpy, these heteroclinic orbits do not form heteroclinic loops; (c) at 
the faces X \J Z where the orbits are heteroclinic orbits connecting singular points located at the 
edge TZxz, these heteroclinic orbits do not form heteroclinic loops. 

edge TZpy. Moreover, given two orbits 71 C 3^ and 72 C S such that a;( 7 i) = 0 ( 72 ), 
then U)('y 2 ) / «(7i)- 

fcj Each of the two limit sets of every orbit in XL)Z is a singular point contained in the 
edge TZxz- Moreover, given two orbits 71 C Af and 72 C .Z such that a;( 71 ) = 0 ( 72 ), 
then U)('y 2 ) / «(7i)- 

Therefore when k G VS-f. (7 iS'*' or k G VS- 7 S~, the behaviour of the trajectories 
in the whole region T is represented in Figure [3j Moreover since a change in the sign of 
the parameter k has the same effect than a change in the sign of the time variable, see 
system m, the behaviour of the trajectories when k G VS- 7 5'*’ or k G VS^ 7 S , 
follows by changing the direction of the flow in Figure [3j 

From Theorem [L2] we conclude that the bifurcation taking place at the manifold S is 
not only characterized by the behaviour of the flow in the interior of T. In addition, it 
must be described by taking into account the changes of the limit sets Sxz and Spy at the 
boundary of T. Hence when k G 7^5+ 7 iS"*" the orbits in the faces X U Z are organized 
in spirals around the segment Sxz moving away from it; and the orbits in the faces T U S 
are organized in spirals around the segment Spy approaching it. When k G VS 7 S, the 
segment Sxz reduces to the singular point ipxz and the segment Spy reduce to the singular 
point furthermore the flow in the faces XCZ and TUS describes heteroclinic orbits 
around them. Finally, when k G P5+ 7 S~ the orbits m. X \J Z are organized in spiral 
around the segment Sxz approaching it; and the orbits in the faces TUS are organized in 
spirals around the segment Spy moving away from it. From this we denote the bifurcation 
taking place at the manifold VSCS by a focus-center-focus bifurcation. The bifurcation 
set of system ([2]) is drawn in Figure 01 

The paper is organized as follows. In Section 2 we analyze the existence and the local 
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Figure 4. Representation of the bifurcation set in a two dimensional parameter space. 

behaviour of the singular points both in the interior and in the boundary of T. In Section 
3 we deal with the first integrals of the flow and we characterize the integrability of the 
flow. Using these first integrals, in Section 4 we analyze the flow at the boundary of T. 
In Section 5 and by using again the first integrals we analyze the flow in the interior of 
T and we prove the main results of the paper. 


2. Singular points 


In the following proposition we summarise the results about the existence, location and 
stability of the singular points of system Q. 

Proposition 2.1 The half straight lines TZpy and TZxz o,re formed by singular points. 

(a) //k G J\fZ there are no other singular points in the boundary of the simplex. 

(a-1) Suppose that k € VS n<S. The open segment 


R = 


ks ki - {ki + k^)z 
ki ' k4 + ki ’ 


0 < z < 


ki 


ks + k4 


is formed by all the singular points in the interior of the region T. Moreover 
the Jacobian matrix of the vector field evaluated at each of these points has one 
real eigenvalue equal to zero and two purely imaginary eigenvalues. 

(a-2) Suppose that k G MZ \ {VS fl S} . There are no singular points in the interior 
of region T. 

(b) Suppose that k 0 MZ and k / 0, in this case there are no singular points in the 
interior ofV. In fact the singular points are on the boundary ofT and they complete 
either edges or whole faces. 
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Proof: Straightforward computations show that the half straight lines TZpy and TZxz are 
formed by singular points. 

(a-1) Suppose now that k G MZ. Hence none of the components of the parameter k is 
zero. In this case the singular points are given by the solutions to the following systems 


X = 0 
yz = 0 > 

z{-k 2 y + kail -y- z)) = 0 ^ 


—x(l — X — z) = 0 
y = 0> 

z(l — X — z) = 0 


x(kiy - k 4 (l -x-y)) = 0 
—yx = 0 > 
z = 0 

/ 


kiX + {ki + kify + kiZ = k^ 
{ki + k4)y + {ka + k^fz = k^ > 
{k2k4 - kika)iy + z) = ^ 2^4 - kika ^ 


where in the last one we impose xyz 7 ^ 0 to avoid repetitions. From the three first systems 
it is easy to conclude that there are no other singular points than those in the half straight 
lines TZpy and TZxz- With respect to the last one we distinguish two situations. 

First let us suppose that k 0 5, that is ^ 2^4 — kika 7 ^ 0. From the third equation it 
follows that y + z = 1, and therefore a; = 0. Since A:iA :4 7 ^ 0 from the first equation we 
conclude that y = 0 and z = 1. Hence the singular point is one of the endpoints of the 
edge TZxz', i-e. it does not belong to the interior of T. 

Suppose now that k G 5, that is k 2 k 4 — kika = 0. Thus the linear system is equivalent 
to the following one 


k 4 X + (ki + k/i)y + k/iz = /C 4 
(/ci + kjfjy + ika + kjfjz = k4. 


If ki + k 4 = 0, then from the first equation we obtain x + z = 1. Therefore y = 0 
and the singular point belongs to TZpy. On the contrary, if ki + k 4 7 ^ 0, then there 
exists a straight line of singular points parametrically defined by x = zka/k 4 and y = 
{k 4 — (ka + k 4 )z)/iki + k 4 ). Since the singular points in the interior of T must satisfy 
that x> 0 , y> 0 ,z >0 and x + y + z < 1, then there exist singular points in the interior 
of T if and only if 

^ > 0 ^3 + ^4 ^ ^ k4 ^ / ka + k4 \ ^ fci ^ ^ ^ 

/c4 ’ /ci + k4 ki + k4^ k4 \ki + k4 J ki + k4^ 


It is easy to check that the previous inequalities are equivalent to 


^>n ^3 + k4 k4 

k4 ' ki + k4 /ci + A:4 ’ 



z > 0 . 


Since k G 5 we have ki/k 4 = ^ 2 /^ 3 - Therefore we conclude that there exist singular 
points in the interior of T if and only if all the components of k have the same sign; that 
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is k G DS. In such case these singular points are given by 

ks k4 - {ks + k4)z r, ^ ^ ^4 

^ = t. + fa • 


( 4 ) 


which proves statement (a-1). 

(a-2) The Jacobian matrix of the vector field defined by the differential equation ([2]) 
and evaluated at the singular points ([!]) is given by 


( ksz {k2 + k3)z ksz ^ 

-kiy 0 k 2 y 

-ksz -{k2 + k3)z -k3Z 

The characteristic polynomial is equal to A(A^ + 6) = 0, where b = zy{ki + k 2 ){k 2 + fcs). 
Since k G VS the coefficient b is positive. Then we get one zero eigenvalue and a pair of 
complex conjugated eigenvalues with zero real part. 

(b) If k 0 J\fZ and k 0, then at least one of the coordinates of k is equal to zero and 
at least one is different from zero. Without loss of generality we suppose that /ci = 0 and 
^2 7^ 0. From the second equation in ([2|) it follows that the coordinates of the singular 
points satisfy yz = 0. Therefore the singular points are contained in the boundary of T. 
Moreover from the remainder equations in (l2|) it follows that —k 4 x{l — x — y — z) = 9 
and k 3 z{l — x — y — z) = 9. We conclude that, depending on whether the parameters fes 
and ki are zero or not, singular points complete either whole faces or edges, respectively. 


In the next result we deal with the singular points located at the edges TZpy and TZxz 
which are not on the segments Spy and Sxz, respectively. Note that these points are not 
hyperbolic singular points, so that we can not apply Hartman-Grobman Theorem to 
describe the behaviour of the flow in a neighbourhood of them. 

Proposition 2.2 // k G VS \ S, then no singular point in TZpy \ Spy and TZxz \ Sxz is 
the limit set of an orbit in the interior of the region V- 

Proof: Let p be a point in the set TZpy \ Spy, that is p = (xq; 0,1 — xq) where either 


xq > max 


k2 _ ks \ 

fci + A:2 ’ ^3 + ^4 i 


or xq < min 


k2 k3 \ 

fei + A:2 ’ fca + /c4 J ’ 


( 5 ) 


see expression ([3|). If we consider a point p in the set TZxz \ Sxz, the following arguments 
can be applied in a similar way. 

Through the change of variables x = x — xq, y = y and z = z — 1 + xq, system m can 
be written as system x = + Q(x) where x = (x, y, z)"^, 


A = 


^ kiXo (fci + A:4 )xo A:4Xo 

0 k2- {ki + k2)xo 0 

\k3{xo - 1) {k2 + A:3)(xo - 1) k3{xo - 1) 
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and 


Q(x) 


^ x{kiX + {ki + ki)y + k^z) ^ 
y{k 2 Z - kix) 

\z{-k3X - {k2 + k3)y - k^z) 


The eigenvalues of the matrix A are Ai = 0, A 2 = {k 3 +k 4 )xo—k 3 and A 3 = ^ 2 —a;o(A:i+A: 2 ). 
From ([5]) it is easy to conclude that A 2 A 3 < 0. Therefore there exists a regular matrix P 
such that PAP~^ = diag{h, A 2 , A 3 }. 

Going through the change of coordinates Xp = Px the system can be rewritten as 


{kihxp + ki{xo - l){k3 + k4)yp + k4{xo - l){k2 + k4)zp) 
Vp = ((-^2 - {k3 + k4)xp)kik4Xo + - xq) + klxo)yp 


< +k4{kik4Xo + k2k3{l - xo))zp) 


( 6 ) 


^p = ((-^3 + {k2 + ki)xp)kik4Xo - ki{kik4XQ + ^ 2 ^ 3(1 - xo))yp 

-k2ik'^{l - xq) + klxo)zp) 

System ([6]) has two invariant planes {yp = 0} and {zp = 0} intersecting at a straight line 
formed by singular points, which corresponds to the segment TZpy. The direction of the 
vector field in a sufficiently small neighbourhood of the origin satisfies that 

sign(yp) = sign(yp)sign(A2) 
sign(ip) = sign(zp)sign(A3). 

We conclude that the origin is neither the a-limit set nor the ca-limit set of any orbit 
in the interior of the regions {yp > 0,Zp > 0},{yp > 0,Zp < 0},{yp < 0, Zp > 0} and 
{yp < 0, Zp < 0}. From this we conclude the proposition. ■ 


3. Invariant algebraic surfaces and first integrals 


In 1878 Darboux showed how to construct first integrals of a planar polynomial vector 
field possessing sufficient invariant algebraic curves. Recent works improved the Dar- 
boux’s exposition taking into account other dynamical objects like exponential factors 
and independent singular points, see 0], and for more details. The extension of 
the Darboux theory to n-dimensional systems of polynomial differential equations can 
be found in the work by Llibre and Rodriguez 1^. A brief introduction to the three 
dimensional case can be found in 
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Following a first integral of system ([2]) is a real function F non-constant over the 
region T and such that the level surfaces Fc = {{x,y,z) € T : F{x,y,z) = C} are 
invariants by the flow; that is 


OF. OF. dF . ^ 

XF = —X + —y + — z = 0 
ox oy oz 


where X = x-^ + y-^ + z-^ is the vector field associated to the system of differential 
equations. Thus the existence of a first integral allows the reduction of the dimension of 
the problem by one. Moreover, the existence of two independent first integrals allows the 
integrability of the flow. 

Let / € M[x, y, z] be a polynomial function. The algebraic surface / = 0 is called an 
invariant algebraic surface of the system ([2|) if there exists a polynomial K G R[a:, y, z] 
such that Xf = K f. The polynomial K is called the cofactor of /. The following result 
is a corollary of Theorem 2 in Q] . 

Theorem 3.1 Suppose that the polynomial vector field ^ admits p invariant algebraic 
surfaces fi = 0 with cofactors Ki for i = 1, 2,... ,p. If there exist € M not all zero 
such that KKi = 0, then the function fi^f^^ ■ ■ ■ fp” is a first integral of the vector 
field 

Now we deal with the existence of Darboux type first integrals of system ([2]). Consider 
the algebraic surfaces /i(x, y, z) = x, f 2 {x, y, z) = y, f^{x, y, z) = z and f^ix, y,z) = x + 
y + z — 1. It is easy to check that Xfi = fiKi, with i = 1, 2, 3,4, where iLi(x, y, z) = k^xF 
{ki + ki)y + k^z-ki, K2ix,y,z) = -kix + k2Z, K^{x,y,z) = -k2,x-{k2 + k^)y-k2,z + k^ 
and Ki{x,y,z) = k^x — k^z. Therefore /i = 0 is an invariant surface with cofactor Ki, 
with i = 1 , 2 ,3,4. 

From Theorem 13.11 if there exist Aj not all zero and such that Ylt=i = 0, then 
F = is a hrst integral of system (l 2 |). Since 


4 

\Ki = {X^ki - \2ki)x + {Xiki - \2,k2)y 

i=l 

+ {X 2 k 2 - Xik 2 ,)z + (Aa/sa - Ai/c 4 )(l - x - y - z) 


the existence of such Aj is equivalent to the existence of non-trivial solutions of the 
homogeneous linear systems 


/ ki -k2\ 

[-k4 ks ) [XsJ 



( k2 -ks\ f A2\ 
[-k, k4 J 1^A4 J 



( 7 ) 


Note that the determinant of both previous systems is equal to kik^ — ^ 2 ^: 4 . Therefore 
when k belongs to the set <S there exist Darboux type first integrals of system ([2|) . 

Under the assumption k G 5 the linear system ([7|) has the following non-trivial solu¬ 
tions (Ai, A 2 , Aa, A 4 ) — >■ (/c 2 ,0,/ci, 0), (0, A: 3 ,0,/C 2 ), (A: 3 ,0,/i; 4 ,0) and (0, A: 4 ,0, fei). Therefore 
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the functions H{x,y,z) = V{x,y,z) = y^^{l — x — y — z)^'^, H{x,y,z) = x^^z^'^ 

and V = y^‘^(l — x — y — z)^^ are first integrals. In fact 


XH = x^'^z^^{l — X — y — z){kik'i — A; 2 ^ 4 ) 

XV = y^^(l -x-y- z)^^x(k 2 k 4 - kiks), 

XH = x’^^z^^ikiks - k2k4)y ^ ^ 

XV = y^^{l — X — y — z)^^ (A; 2^4 — kik^)z 

which vanish in the whole region T only when k G <S. 

Proposition 3.2 Consider the functions H{x,y,z) = x^'^z’^V H{x,y,z) = x^^z’^*, 
V (x, y, z) = y^^ {1 — X — y — z)^^ and V(x, y, z) = y^^{l — x — y — z)^^. 

(a) IfV€S r\J\fZ, then H, V, H and V are first integrals which satisfy that 
and V^^ = V^'^. Moreover H and V are independent. 

(b) //k G iS \J\fZ, then two of the previous functions are first integrals and they are 
independent. 

(c) // k 0 5, then none of the previous funetions is a first integral in T. 

Proof: (a) Consider that k G 5 fl J\fZ. Since every coordinate of k is different from 
zero it follows that H, V, H and V are not constant in T. Therefore all of these func¬ 
tions are first integrals. It is easy to check that and V^^ = V^'^. More¬ 
over since XH{x,y,z) = {k2Z,0,kix) and XV{x,y,z) = — x — 

y — (—fc 2 y, ^ 3(1 — X — y — z) — k2y^ —k2y ), both integrals are dependent only on 

points satisfying ^ 3(1 — x — y — z) = k 2 y and k 2 Z = kix. Taking into account that ^2 7 ^ 0 
it follows that this set has zero Lebesgue measure. Then H and V are two independent 
first integrals. 

(b) Consider now that k G S\MZ. Hence k has one coordinate which is different from 

zero. Without loss of generality we assume that ki 7 ^ 0, the remainder cases follows in 
a similar way. It is easy to check that H and V are not constant in T, and therefore 
they are first integrals. Since VH(x,y, z) = (k2Z,0, kix) and VV{x,y, z) = 

y(^'‘-i)(l — x — y — {—kiy, ^ 4(1 — x — y — z) — kiy, —kiy ), both integrals may be 

dependent only on points satisfying /c 4 (l — x — y — z) = kiy and k 2 Z = kix. Therefore H 
and V are independent. 

(c) The statement follows straightforward form expression ([ 8 ]). ■ 


4. Behaviour at the boundary 

As we have proved in Proposition 13.21 some of the functions H, H, V and V are first 
integrals over the whole region T only when k G <S. Nevertheless the restriction of these 
functions to a particular face of T results in a first integral even when k 0 5. In fact, 
denoting by H\y the restriction of the function H to the face T, from expression Q it 
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Figure 5 . (a) Foliation over the face y defined 
C2 < C3 < C* and k in VS \ S. (b) Foliation 
or 14- = or % = where 0 < 
figure (b) is represented in {u, u)-coordinates. 



by the level curves H\y = where 0 < Ci < 
over the corresponding face of the level curves 
Cl < C2 < C3 < C* and k in VS \ S. Note that 


follows that XH\y = 0. Therefore the level curves H\y = are invariant by the flow. 

Under the assumption > 0, these level curves define a foliation of y whose leaves 

are given by the arcs of hyperbolas \ z = Cx *“ 4 1 where 

1 J 0<C<C* 


kj ( ks \ '=4 
k4 ks \kii -\- kj,) 


(9) 


Furthermore, every leaf with 0 < C < C* intersects the segment TZpy at exactly two 

points, see Figure [5ja). The value C = C* leads to a unique intersection point with 

coordinates x = k 3 /{k 3 + ^ 4 ) and 2: = k ^/(fca + k^). Since in the face y we have y = 0, it 

follows that the point corresponding to C* is the point (\py defined in ([3]). 

Similarly, the restriction of V, V and H to the faces T, Z and S respectively, are 

first integrals even when k ^ 5, see expression ([8|). Consider the changes of variables 

{u,v,a,f3) —>■ {y,z,k 2 ,k 3 ),{y,x,—ki,—ki) or {x,y,ki,k 2 ), depending on the face X, Z 

or S we are looking at. Under the assumption af3 > 0, the level curves U|.y = C^^Vlz = 

and define a foliation on the corresponding face, whose leaves are given 

by the unimodal curves \ v = 1 — u — Cu~^ \ where 

I i 0<C<C* 

C*- ^ 

q: + /3 \Oi P J 

Every leaf with 0 < C < C* intersects the segment {u = 0,0 < u < 1} at exactly 
two points, see Figure [Hb). The value C = C* leads to a unique intersection point 
{f3/{a + f3),0) . Going back through the change of variables and adding the variable 
which does not appear in such change, that intersection point coincides with qj,^, ^pxz or 
Ppy depending on the change of variables. 

Using the geometric information of the aforementioned foliation, in the next result we 













March 20, 2015 0:18 Dynamical Systems 3DLV'DSIJ2 


14 


Taylor & Francis and I. T. Consultant 


summarise the behaviour of the flow of system ([ 2 ]) at the boundary dT for k € VS. 
Lemma 4.1 

(a) //k G VS, then each of the two limit sets of every orbit contained in 3 ^U S (respec¬ 
tively, XC Z) is formed by a singular point contained in the edge TZpy (respectively, 
in the edge TZxz)- 

(b) //k G VSriS, then for every pair of orbits 71 C 3 ^ and 72 C S (respectively, 71 C if 
and 72 C .ZJ satisfying that a;(7i) = 0(72), it follows that 0(71) = a;(72). 

(c) //k G VS\S, then for every pair of orbits 71 C 3 ^ and 72 C S (respectively, 71 C if 
and 72 C ZJ satisfying that a;(7i) = 0(72), it follows that 0(71) 7 ^ ^(72). 

Proof: We restrict ourselves to consider orbits in the faces 3iUS. The study of the orbits 
in the faces X U Z follows in a similar way. 

(a) Suppose that k G VS. Hence ^3/04 > 0. Therefore every orbit 71 in y is contained 
in a leaf of the foliation z = Cx with 0 < C < C*, which is an arc of a hyperbola 
intersecting the edge TZpy at exactly two points. Since there are not other singular points 
in y, see Proposition I2.11 ai. we conclude that each of the two limit sets of 71 is one of 
these intersection points. 

On the other hand we have ^ 2^1 > 0. Therefore every orbit 72 in S is contained in a 
leaf of the foliation y = 1 — x — Cx *^ 1 , which is an unimodal curve intersecting the edge 
TZpy at exactly two points. We conclude again that each of the limit sets of 72 is one of 
these intersection points. 

(b, c) Taking into account that S is given by the relation z = l — x — y, we express 
the leaves in S as a function z{x) in the following way z = Cx ‘’i. 

Let p = (xo,0,1 — xo) be a point in the edge Ppy. There exist two positive values Ci 

fcs ^2 

and C 2 such that both the leaf z = Cix *=4 in the face y and the leaf z = C 2 X '=1 in 
the face S contain the point p. On the other hand the leaf in the face y intersects TZpy 
at a new point (xi,0,1 — xi) and the leaf in the face S intersects TZpy at a new point 
{x 2 , 0,1 — X 2 ). Since two hyperbolas either intersect at most at one point or they coincide, 
we conclude that k^^ki = k/^k 2 if and only if xi = X 2 . ■ 


5. Behaviour in the interior 

In this last section we deal with the proof of the main theorems of the paper. The next 
result is a technical lemma which describes planar flows under integrability conditions. 

Lemma 5.1 Let x = f (x) he a planar system of differential equations and let U be a flow- 
invariant region in Assume that both the boundary of U is formed by a heteroclinic 
loop and there exists exactly one singular point p contained in the interior of U. If there 
exists a first integral H defined over U which is non-constant over open sets, then every 
orbit in the interior of U but the singular point is a periodic orbit. 
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Proof: Let U denote the interior of the region U. It is easy to check that any orbit 7 in 
U \ {p} has its limit sets contained in the boundary of U \ {p}. Otherwise first integral 
H would be constant over one of the open regions limited by 7 or over the whole open 
region U \ {p}. Hence there are not homoclinic orbits to the singular point p. From the 
Poincare-Bendixson Theorem, at least one limit set of 7 is a periodic orbit Pi surrounding 
p and contained in f7 \ {p}. 

Applying now similar arguments to the other limit set, it follows that this limit set is 
also a periodic orbit r 2 contained mU \ {p} and surrounding p. Since H is not constant 
over open sets we conclude that Fi and r 2 are the same periodic orbit and this periodic 
orbit coincides with 7 . Therefore every orbit in IJ \ {p} is a periodic orbit. ■ 

Proof of Theorem HTJ (a) Under the assumption k G VS PI S, system ([2]) is 
integrable and the functions H and V are two independent first integrals, see Proposition 
l3.2lf aL Since any level surface Vc is invariant by the flow, we can consider the restriction 
of the flow to each of these surfaces. Of course this restricted flow is also integrable 
because the restriction of the function V to Vc is a first integral. On the other hand 
there exists exactly one singular point in the interior of T-Lci which comes from the 
intersection of the manifold T-Lc ^-Rcl the segment R defined in the Proposition 12.ll fa-ll. 

The map 7r{x, y,z) = {x,y) projects the manifold PLc over a compact region U. More¬ 
over the Jacobian matrix of tt at the point 7 r“^(x,y) defines a flow over U which is 
differentially conjugate to the restricted flow over T-Lc- Since there exists exactly one 
singular point in the interior of U, the function V ott~^ is a first integral over U which is 
non-constant over open sets and the boundary of U is formed by a heteroclinic loop (see 
Lemma l4.1f bll. from Lemma l5.II it follows that every orbit in U but the singular point 
is a periodic orbit. Therefore, every orbit over Re but the singular point is a periodic 
orbit. This result is independent on the level surface we are working at, hence every orbit 
in the interior of the region T, but the singular points, is a periodic orbit. 

The behaviour of the flow at the boundary of T when k G VS H S can be obtained 
from Lemma iJ.H bl. 

(b) Consider now that k 0 VS H S and k 7 ^ 0. We distinguish between two situations: 
first we suppose that k G 5 \ VS. In such case k belongs to the manifold S. From 
Proposition 13.21 it follows that at least one of the functions H, U, LJ or U is a first integral. 
Without loss of generality we can assume that iL is a first integral. Hence any level 
surface Rc is invariant by the flow and we can consider the restriction of the flow to 
Rc- From Proposition 12.11 there are not singular points in the interior of Rc- Applying 
the Poincare-Bendixson Theorem to the flow in the level surface Rc, we conclude that 
the flow goes from the boundary of Rc to the boundary of Rc- Since these arguments 
are independent on the level surface, it follows that the limit sets of every orbit in the 
interior of T is contained in dT- 

Suppose now that k 0 5 and k 7 ^ 0. Since one of the coordinates of k is different from 
zero, the level surfaces of at least one of the functions H, V, H and V can be expressed 
as the graph of an explicit differentiable function. For instance if 7 ^ 0, then Rc^^i is 




March 20, 2015 0:18 Dynamical Systems 3DLV'DSIJ2 


16 


Taylor & Francis and I. T. Consultant 


z z 


Vc 


Figure 6. Representation when k € 7 ^<S+ HiS''' of: the positive invariant regions (in grey) limited 
by the level surfaces He and He] the negative invariant regions (in grey) limited by the level 
surfaces Ve and Vc; the segment Spy formed by the w-limit set of each orbit in the interior of F ; 
and the segment Sxz formed by the a-limit set of each orbit in the interior of T. 

_ia. 

the graph of the function z = Cx defined over the face Z. Each of these level surfaces 
split the interior of T into two disjoint connected components. On the other hand since 
k ^ 5 these level surfaces are not invariant by the flow, see Proposition I3.2r c). In fact 
the flow is transversal to them and the direction of the flow through them depends on 
k € iS"'' or k G see expression ([HI). Since as C tends to 0 or to C* the level surfaces 
Hc^i tend to the bonndary of T, we conclude that the flow in the interior of T goes 
from one part of the boundary to another part of the boundary. That is the limit sets of 
every orbit in the interior of T are contained in dT. 

Thus, in both cases the limit sets of every orbit in the interior of T are contained in 
dT- Since there are no isolated singular points in dT, see Proposition 12.11 we conclude 
that these limit sets are not periodic orbits. ■ 




Note that in the previous proof we have only used that trajectories cross the level 
surfaces of some of the functions H, V, H or V, always in the same direction. This 
suffices to conclude that the limit sets of the orbits in T are contained in the boundary. 
To prove Theorem 11.21 we need to be more precise in the location of these limit sets. To 
reach this goal we will control the geometry of the level snrfaces. 

Proof of Theorem 11.21 (a) Suppose that k € TS+ nS~^. Since k 0 5 the functions 
H and H are not first integrals and each level surface He and He splits the region T 
into two disjoint regions in such a way that the flow goes from one to the other. In fact 
since > 0 the intersection of He with any plane {y = yo}, where 0 < yo < 1) is 
an arc of a hyperbola in the (x, 2 ;)-plane, see Figure 0 Similarly, since kik 2 > 0 the 
intersection of He with any plane {y = yo}) where 0 < yo < 1, is an arc of a hyperbola 
in the (x, zj-plane. The flow through He and through He has the same orientation as 
the vectors VH and VH respectively, see expression ([8]). Since fc, > 0 the coordinates of 
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the gradient VT/ are non-negatives. Hence Vi7 is oriented towards the region containing 
the point ppy, see the shadowed region in Figure [H In a similar way, the gradient V77 
is oriented towards the region containing the point c[py, see Figure [6l Therefore the flow 
evolves from the region containing the origin to the region containing the segment Spy, 
see Figure 0 On the other hand points in TZpy \ Spy are not limit set of orbits in the 
interior of T, see Proposition 12.21 We conclude that the w-limit set of any given orbit in 
the interior of T is a singular point contained in the segment Spy. 

Since k 0 5 the functions V and V are not first integrals. Moreover each level surface 
Vc and Vc splits the region T into two disjoint regions in such a way that the flow goes 
from one to the other. The flow through these surfaces has opposite direction to that of 
the gradients W and VF, see expression ([8]). Since ki > 0 the gradient VF is oriented 
towards the region containing the point and the gradient VF is oriented towards 
the region containing the point Pxz- Therefore the flow in the interior of T evolves from 
the shadowed region in Figure [6] towards the region containing the point (0,0,1). Since 
points in 7Zxz \ Sxz are not limit set of orbits in the interior of T, see Proposition 12.21 we 
conclude that the a-limit set of any given orbit in the interior of T is a singular point 
contained in the segment Sxz-, see Figure 0 

As we have just proved when k G VSj^ ntS"*" the w-limit set and the a-limit set of any 
given orbit in the interior of T is contained in the segments Spyy and Sxz^ respectively. 
Similar arguments apply when k G VS- n S~. 

Note that a change of the sign of the parameter k is equivalent to a change in the sign 
of time in the system of differential equations ([2|). Therefore the behaviour of the flow 
in cases k G 'P5+ H S~ and k G VS- H follows from the cases described above by 
changing the orientation of the orbits. 

(b,c) The behaviour of the flow at the boundary can be obtained from Lemma Id.H cl. 
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